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ABSTRACT 

In magnetic cataclysmic variables (CVs), X-ray emission regions are located 
close to the white dwarf surface, which is expected to reflect a signihcant fraction 
of intrinsic X-rays above 10 keV, producing a Compton reflection hump. However, 
up to now, a secure detection of this effect in magnetic CVs has largely proved 
elusive because of the limited sensitivity of non-imaging X-ray detectors. Here 
we report our analysis of joint NuSTAR/XMM-Newton observations of three 
magnetic CVs, V709 Cas, NY Lup, and V1223 Sgr. The improved hard X- 
ray sensitivity of the imaging NuSTAR data has resulted in the first robust 
detection of Compton hump in all three objects, with amplitudes of or greater 
in NY Lup, and likely <1.0 in the other two. We also confirm earlier report of a 
strong spin modulation above 10 keV in V709 Cas, and report the hrst detection 
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of small spin amplitudes in the others. We interpret this as due to different 
height of the X-ray emitting region among these objects. A height of ~0.2 white 
dwarf radii provides a plausible explanation for the low reflection amplitude of 
V709 Cas. Since emission regions above both poles are visible at certain spin 
phases, this can also explain the strong hard X-ray spin modulation. A shock 
height of ~0.05 white dwarf radii can explain our results on V1223 Sgr, while 
the shock height in NY Lup appears negligible. 

Subject headings: Novae, cataclysmic variables — X-rays: binaries 


Introduction 


In magnetic cataclysmic variables (CVs), a strongly magnetic (~l-200 MG) white dwarf 
primary accretes matter from a Roche-lobe hlling secondary. Intermediate polars (IPs) are 
a subset of magnetic CVs in which the spin period is signihcantly shorter than the orbital 
period. When the magnetically controlled accretion flow hits the white dwarf surface, a 
strong standing shock forms. The predominant cooling mechanism of shock-h eated plasma 


fr K 

in IPs is the collisionally excited, thermal, multi-temperature X-ray emission flAizu 1 119731) . 
Since the shock is formed near the stellar surface, we expect that roughly half the intrinsic 
emission is directed towards the primary and is reprocessed and/or reflected. The strongest 
spectral signatures of reflection are the 10-30 keV Compton reflection hump and the Fe K 
fluorescent line at 6.4 keV. While the latter are routinely detected from magnetic CVs, they 
can also originate from intrinsic absorbers. The importance of reflection has nevertheless 
been inferred thro ugh a statistical analy sis of the equivalent width of the 6.4 keV line vs. 
absorbing column flEzuka fc Ishida II1999I) . 


A direct detection of the Compton reflection hump in magnetic CVs has largely been 
elusive to date, due to a lack of X-ray s pectra of sufficient qua l ity above 10 keV. Posi¬ 


tive detections ar e reported for AM Her flRothschi 


Matt et ah 112000 ) and for EF Eri ( Done et ah 


Done et ah 


d et ah 1981: Beardmore et ah 1995 


19951). How ever, the studies on AM Her used 


(119951) also started with this model, and 


single-temperature models as the basis; 
although they also explored reflection of the expected multi-temperature emiss ion, the data 


were not good enough to distinguish between the two. iBeardmore et ah (120001) found resid¬ 


uals around the Fe K complex in their study of V1223 Sgr using Ginga data, and interpreted 
it as due to reflection. However, they were unable to verify this through the detection of 
the characteristic spectral curvature of reflection on the higher energy side of the 10-30 keV 
hump. 

































3 


Previous hard X-ray data on these magnetic CVs generally had poor signal-to-noise 
ratios and suffered from systematic uncertainties of background models. With the advent of 
NuSTAR, the hrst mission to extend the energy range for imaging X-ray observations up to 
~79 keV, we can obtain much improved hard X-ray spectra. It also enables a reliable study 
of spin modulation of X-rays above 10 keV. In thi s Letter, we present th e initial results from 


our joint observatio ns of 3 IPs using NuSTAR (iHarrison et ah 1201311 and XMM-Newton 


fjjansen et ah Il200l|l . and report the unambiguous detection of Compton reflection hump in 
all three targets. 


2. Observations and Data Reduction 


We have selected V709 Cas, NY Lup, and V1223 Sgr as our targets fo r this study, because 


they h ave the highest observed 14-195 keV Swift/BAT fluxes according to lBrunschweiger et ah 


(120091) . We have observed them jointly with NuSTAR and XMM-Newton during the summer 
of 2014, always with periods of simultaneity; the observations are summarized in Tabled! In 
this Letter, we present only the NuSTAR and the XMM-Newton EPIC data. 


NuSTAR consists of two co-aligned telescopes and the two focal plane modules, FPMA 
and FPMB. We reduced the NuSTAR data using the NuSTAR Data Analysis Software 
(NSuTARDAS vl.4.1) available as a part of HEASOFT v6.16 and latest CALDB hies. The 
NuSTAR observations were taken during intervals of normal Solar activity, and we used 
standard hltering to remove periods of high background during South Atlantic Anomaly 
(SAA) passages and Earth occultation. For each source, we used a 70" radius circular region 
around the known position, and a 100" radius, circular, source-free background region on the 
same detector. We generated the spectra, corresponding response hies and light curves using 
nuproducts. We grouped all spectra using grppha to have at least 25 counts per bin, and 
used the statistic in spectral htting. We produced background-subtracted, barycenter- 
corrected light curves in the 3-10 keV and 10-30 keV energy bands for all three IPs. 


We processed the XMM-Newton data with SAS version 13.0.0 using the latest cali¬ 
bration hies available in November 2014. All observations were performed with the EPIC 
cameras (PN and two MOS) in small window imaging mode to avoid pile-up, with the thin 
hlter applied (with the exception of V709 Cas, where the medium hlter was used). We con- 
hrmed that none of the observations was ahected by pile-up following the thread procedur^. 
Then, we extracted the source spectrum from a circular region of radius 40", centered at 
the known source position. The background was obtained from a nearby source-free region 


^ http://xmm.esac.esa.int/sas/current/documentation/threads/epatplot.shtml 
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of the same size for PN and 150" for the MOSs. For the spectral analysis, we cleaned the 
observation of V709 Cas that was affected by backgronnd flares, by applying an intensity 
threshold, thereby retaining only the hrst ~ 22 ks. We gronped spectral channels with the 
specgroup tool so as to have a minimnm of 25 and 15 connts per bin for the PN and MOSs, 
respectively. For the timing analysis, we applied the barycentric correction nsing barycen 
and produced background subtracted light curves in the 0.3-3, and 3-10 keV ranges using 
the task epiclccorr. 


3. Results 


3.1. Phenomenological Fits 


Initially, we concentrated on what we could infer purely from the NuSTAR data. For 
this, we excluded the Fe K region (5.5-8.0 keV), and the data below 4.5 keV, and ap¬ 
proximated the interstellar plus intrinsic absorber using a single phabs model. To £t the 
continuum, we used 3 different models with different degrees of curvature: a power law, a 
cooling flow model (mkcflow), and a single-temper ature Bremsstrahlung. We the n modeled 
reflection using the convolution model ref lect flMagdziarz fc Zdziarski Il995l) . with the 


abundances hxed to solar as determined by Anders fc Grevesse f 19891) . We also hxed the 
inclination angle of the reflecting surface at the default value of cos /i = 0.45 (/r ~ 63.3°) 
(see §3.2). We htted the data once with the reflection amplitude hxed to 0.0, and once with 
the amplitude thawed. 


In Table[2l we report the best-fit xl and the degree-of-freedom of the latter case, the 
difference in (Ay^) compared to the no rehection case, and the best £t rehection amplitude 
with the 90% conhdence error range. The choice of continuum model affects the Ay^ values, 
with the Bremsstrahlung continuum presents the stringent test of rehection (the smallest 
Ay^). Nevertheless, rehection is detected with high signihcance for all 3 models and for all 
3 objects (at a signihcance of 1—6 x 10“^°, for V709 Cas htted with Bremsstrahlung model, 
or higher, using the F-test). We also illustrate the improvement visually in (Figured]) using 
the single-temperature Bremsstrahlung case. The residuals against the rehection-less model 
(middle panels) clearly show the continuum rehection hump, peaking around 20 keV. We 
are not aware of any explanation of the high energy end (~30 keV) of the observed hump 
other than Compton rehection. 


We conclude that, with NuSTAR data alone, we have achieved an unambiguous detec¬ 
tion of rehection in all three IPs. However, intrinsic absorbers in IPs are complex and often 
include components with Nh ~a few times 10^^ cm“^. Such absorbers ahect the observed 
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continuum up to ~6 keV and also create an Fe K edge at 7.0 keV flEzuka fc Ishida Ill999l) . 
We must determine the depth of the Fe K edge due to these Compton-thin absorbers before 
we can place quantitative constraints on the reflection amplitude (see §3.2). 

We also performed simple phenomenological hts to the XMM-Newton EPIC spectra 
in the 5-9 keV region, where multiple emission lines are evident. Using a Bremsstrahlung 
continuum modihed by a single phabs absorber, we clearly detected the H-like (7.0 keV), 
He-like (6.7 ke V), and near neutral, fl uorescent (6.4 keV) lines in all 3 objects, which is 
typical for IPs (IHellier fc Mukai 1120041) . The hrst two indicate the presence of kT~10 keV 
and ~5 keV plasmas, respectively. The equivalent width of the 6.4 keV line was found to be 
105T11 eV, 132±12 eV and 90 ±18 eV, respectively, for V709 Cas, NY Lup, and V1223 Sgr. 


3.2. Fits to Combined Data 


Any absorbers capable of producing an Fe K edge will necessarily absorb soft X- 
rays below 2 keV almost completely. On the other hand, many IPs often exhibit ad¬ 
ditional s pectral complexities belo w 1 keV, such as the soft, blackbody-like component 
see, e.g.. iBernardini et ah 120121) and the complex absorber may be ionized (see, e.g.. 


de Martino et ah 1 120081) . In this initial study, we have therefore opted to £t XMM-Newton 


EPIC (pn, MOSl and MOS2) spectra in the 1.5-10 keV range simultaneously with the NuS- 
TAR (both modules) spectra in the 3.5-78 keV range. Our preliminary analysis down to 0.3 
keV conhrms that the choice of low energy threshold does not signihcantly affect the depth 
of the Fe K edge in the model. 

For these joint hts, we used two different combinations of models. In the hrst case 
(“2T”), we approximated the underlying multi-temperatu re thermal ernission, evidenced by 
the p resence of H-like and He-like Fe lines, using two mekal ( Mewe et ah Ill985 : Liedahl et ah 
19951) components of diherent temperatures. A single absorber is clearly inadequate to de¬ 
scribe the data, so we approximated the complex absorber using one fully covering and one 
partial-covering absorbers. We applied the reflect model to the 2 mekal models and added 
a Gaussian to represent the 6.4 keV line of cold iroij^. 

In the second case (“mkchow”), we modeled the und erlying multi-temperature plasma 


emission using the isobaric cooling how model, mkcflow (IMushotzky fc Szvmkowiak Ill988 


^In xspec notation, this model can be written as 

constant*phabs (partcov*phabs) (reflect (mekal+mekal)+gaussian) where the constant allows for cross¬ 
normalization constant, and the convolution model partcov is used to create a partial-covering absorber 
model based on phabs. 
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Mukai et al. 1 120031) . Although we expect post-shock plasma to cool to the photospheric 


temperature of the primary, the data above 1.5 keV have no sensitivity to plasma temper¬ 
atures well below 1 keV. We therefore used a fixed minimum temperature of 0.0808 keV, 
the lowest value allowed by the mkcflow model. The intrinsic absorption in IPs is likely to 
occur in the immediate pre-shoc k flow, which is poorly appr oximated by one or several dis¬ 
crete partial-covering absorbers: iDone fc Maedziarz fjlOOa l pointed this out and provided 
a complex absorber model with an power-law distribution of covering fraction as a function 
of Nh, pwab. We used pwab with minimum Nh set to 1.0x10^® cm“^ (effectively 0.0) as well 
as a phabs to model the interstellar absorptiorj^. 


In both these models, we allowed the overall abundances of the thermal emission to 
vary (linked between the two mekal components for the “2T” model), which is necessary to 
reproduce the strengths of the H-like and He-like Fe lines. Moreover, we assumed that the 
reflection was from the part of the white dwarf surface that is covered by freshly accreted 
material, with the same abundances as the X-ray emitting plasma. We initially fixed the 
angle of reflection at cos /i=0.45 as before. 


We allowed the cross-normalization factor relative to NuSTAR FPMA to vary, to al¬ 
low for sour ce variability and for the current limitatio n in the cross-calibration of these 
instruments fiBurrows et al. I2014J: iMadesn et al. l2015h . The best-ht values were all found 
to be within a few percent of 1.0. However, we found several specihc manifestations of 
cross-calibration limitations at, e.g., 7-10 keV. This is the likely reason why the joint NuS¬ 
TAR/ XMM-Newton hts results in with xt > which might otherwise suggest that these 
hts are unacceptable. We summarize the reflection amplitudes and other selected results in 
Table[2l and show the results for V1223 Sgr in Figure[2l 


The derived reflection amplitudes are lower than in NuSTAR-on\j hts, because our full 
hts indicate the presence of a complex intrinsic absorber up to Nh ~a few times 10^^ cm“^, 
which the NuSTAR-only hts did not account for. Correspondingly, these hts result in smaller 
but still statistically signihcant values of Ay^. We also note that there is good agreement in 
the rehection amplitude between the two models for V709 Cas and for V1223 Sgr, and that 
they are signihcantly less than 1.0: ~0.65 (0.45-0.9) for V1223 Sgr and ~0.35 (0.1-0.6) for 
V709 Cas. For NY Lup, the diherent models disagree but the rehection amplitude appears 
to be >0.8, possibly as high as 1.8. Alternatively, if we allow for diherent viewing geometry, 
we can obtain hts with similar quality with rehection amplitudes ~1. For NY Lup htted 
with mkcflow, this can be achieved with ju ~ 30°, resulting in a rehection amplitude of 0.83- 


^This model can be expressed as constant*pliabs*pwab(reflect*mkcflow+gaussian) in xspec nota¬ 
tion. 
















7 


1.20. For V709 Cas and V1233 Sgr, the data can be reconciled with reflection amplitndes of 
~1.0 only with edge-on > 70°) geometries. Finally, we note t hat we do not confirm th e 


previons claim of an ionized Fe K edge at ~8.0 keV in V709 Cas fide Martino et ah Il200ll) . 


3.3. High Energy Spin Modulation 


Energy-dependent X-ray spin modnlation is a definin g characteristic of IPs, and is gen¬ 


erally nnderstood as dne to variable complex absorber flNorton &: Watson Il989l) . This 


mechanism predicts small or negligible spin amplitndes above 10 keV. The only previons 
claim of spin modnlation above 10 keV is for V709 Cas, which shows >10 keV modnla¬ 
tion in phase with the lower-energy m odnlation in the non-imaging RXTE and Beppo-SAX 
(PDS) data (jde Martino et ah Il200ll) . To confirm this earlier resnlts, and to search for snb- 
tle spin modnlation in the other objects, we folded the XMM-Newton and NuSTAR light 
cnrves in three energy ba nds (Fignre[2]) on the k nown spin periods of 312 .75 s for V709 Cas 
f de Martino et ah 2001). 693.01 s for NY Lnp f de Martino et ah 2006), and 745.63 s for 


V1223 Sgr fjOsborne et ah Ill985l) . nsing the integer part of Modified Jnlian Date of the 


starting time of each observation as the epoch of phase 0.0. 

We fonnd the pnlsed fractions of 10-30 keV spin modnlation to be 17.0±1.7%, 11.6±1.6%, 
and 4.2 ±1.1%, respectively , for y 709 Cas, NY Lnp, and V1223 Sgr, nsing the definition 
in, e.g., iBernardini et a.1. (1201211 . Onr resnlt for V709 Cas is consistent with that of 


de Martino et ah (1200 ll) . In NY Lnp, the 10-30 keV spin modnlation is ronghly sinnsoidal 
and in phase with the modulation below 10 keV; The 10-30 keV modulation of V1223 Sgr 
is the weakest and non-sinusoidal (Figure[3]). 


4. Discussion 


The reflection signal is expected to be the strongest normal to the r eflecting surface, even 
before takin g projection effect into accoun t (IGeorge fc Fabian 1 1199111 . with limited energy 
dependence. iMaedziarz fc Zdziarski (119951) chose cos p = 0.45 as the default for their xspec 
model as having a spectral shape close to the ensemble average. In our fiducial fits, we used 
this default and found the reflection amplitude to be <1.0 in V709 Cas and V1223 Sgr, and 
high (quite possibly >1.0) in NY Lnp. However, the primary surface can cover at most 27r 
steradian of sky. Our data on NY Lnp can be reconciled with this if the viewing geometry 
is more face-on {p ~ 30°). This would be the case if the binary inclination (z) is ~ 30° and 
the co-latitude of the emission region (/3) is small, which we believe is plausible for NY Lnp. 








































If we force the reflection amplitude to be 1.0 in V709 Cas and V1223 Sgr, the spin- 
averaged /i must be > 70°. This would require either a nearly aligned rotator seen edge-on 
or equatorial emission regions seen pole-on. However, our data are also consistent with 
non-negligible shock height. Unfortunately, given the abundance and Nh values we infer, 
we cannot make a strong inference on the ir nportance of reflection based on the equivalent 
widths of the 6.4 keV lines (see Figure 7 of lEzuka fc Ishida lllOOOl) . Therefore, we use the 
spin modulations to break this degeneracy. 

The soft X-rays (<10 keV) are modulated predominantly due to photoelectric absorption 
and are at a minimum when the upper pole (the pole on the hemisphere nearer the Earth) 


points most directly towards us flHellier et ah Ill987|l . In contrast, the /i dependence of 


reflection should produce a hard X-ray maximum at this phase, resulting in anti-phased spin 
modulations between the soft and the hard bands. Since this is contrary to what we observe 
(§3.3), the hard X-ray spin modulations have a different mechanism. 

For NY Lup and V1223 Sgr, the observed hard X-ray spin amplitudes are of order 10% 
or less. The complex absorbers we infer for these objects (Table ip may prod uce a spin 
modulation of such amplitudes via Compton scattering (Figure 6 of iRosen II1992I1 . provided 
that the angle-dependent reflection effects are small (i.e., the range of p must be small). 


Finite shock height is another potential mechanism for spin modulation flMukai Ill999l) . 
If the shock height is negligible, as is often assumed in studies of reflection and spin modula¬ 
tions, the horizon is 90° from zenith and the primary surface covers 27r steradian, implying 
a reflection amplitude of 1.0. Moreover, in a pure dipole geometry, one and only one pole 
is observed at all phases, as the disappearance of one behind the primary rim exactly coin¬ 
cides with the the appearance of the other. Once we relax this assumption, the horizon is 
>90° from zenith, the surface covers < 27r steradian, and there is a range of viewing angles 
where the X-rays from both emission regions are simultaneously observable. This leads to 
an energy-independent spin modulation, since this effect is purely geometrical. 

To explain the hard X-r ay spin modulation o f V709 Cas, already kno wn to be non- 
sinusoidal from ROSAT data (iNorton et ah Ill999h . lde Martino et ah (120011) proposed that 
its shock height is ~0.2 white dwarf radii (R^^; see their Figure 9). Such an X-ray source 
has a horizon angle of 123.5°, both poles are observable at the viewing angles of 66.5°- 
123.5°, and the reflection amplitude is 0.45. This scenario provides an explanation for our 
observation that the reflection amplitude is low, the 6.4 keV line is relatively weak, and the 
spin modulation is even above 10 keV. 


The shock in NY Lup appears much closer to the primary surface. As for V1223 Sgr, 
a characteristic height of 0.05 Ri„rf would result in a reflection amplitude of 0.7, consistent 
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with the data assuming a non-edge-on viewing geometry, and a horizon angle of 107.7°. The 
geometry of V1223 Sgr is poorly known; we believe that there are plausible combinations of 
i and {3 (e.g., ~ 50° and ~ 20°, respectively) such that the lower pole is never visible. We 
can then explian the observed spin modulation as due primarily t the Compton scattering 
mechanism. 


5. Conclusions 

We have unambiguously detected the continuum reflection bump in all three IPs we ob¬ 
served jointly with NuSTAR and XMM-Newton. For two of them, V709 Cas and V1223 Sgr, 
either we view the reflecting surface edge-on or the reflection amplitude is low. Considering 
also the high energy spin modulation we detect in V709 Cas, we argue that a post-shock 
region with a signihcant height above the white dwarf surface can explain all available data. 
We plan to continue our analysis using the XMM-Newton RGS and OM data, using EPIC 
data down to the lowest energy, and performing spin-phase resolved spectroscopy of EPIC 
and NuSTAR data, to investigate if this framework can provide a quantitative description 
of the data. 
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Fig. 1.— Bremsstrahlung fits to NuSTAR spectra of 3 IPs with and without reflection. In 
the top panels, data are shown along with best-£t models with (blue) and without (red) 
reflection. The two lower panels show data to model ratios for hts without (middle) and 
with (bottom) reflection. 
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Channel Energy [keV] 

Fig. 2.— Cooling flow fit to the combined NuSTAR and XMM-Newton spectra of V1223 
Sgr. Only NuSTAR FPMA (black), XMM-Newton EPIC PN (red) and MOSl (blue) data 
are shown for clarity. 
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Fig. 3.— Spin-folded light curves in 3 energy bands. Two cycles are shown. Data used 
are obtained with XMM-Newton EPIC PN (shown in red: 0.3-3 keV and 3-10 keV) and 
NuSTAR FPMA (black: 3-10 keV and 10-30 keV). 


















Table 2. Selected Spectral Fit Results 


Object 

Data 

ModeP 

Xl/^OF 


Amplitude Nh® 

10^^ cm“^ 

PCF/a® 

kT'^ 

keV 

Abund® 

V709 Gas 

NuSTAR 

PL 

1.16/557 

78.9 

2.16 (1.57-2.91) 







mkcflow 

1.02/557 

81.7 

1.71 (1.27-2.24) 







Brems 

1.05/557 

41.7 

0.80 (0.56-1.07) 






Joint 

2T 

1.29/1094 

8.7 

0.27 (0.11-0.44) 

1.5 

31% 

27.0 (28.6) 

0.34 



mkcflow 

1.30/1096 

12.4 

0.38 (0.21-0.60) 

4.7 

-0.84 

46.6 (48.2) 

0.28 

NY Lup 

NuSTAR 

PL 

1.02/401 

219.4 

3.95 (3.07-5.12) 







mkcflow 

0.97/401 

279.4 

2.68 (2.13-3.26) 







Brems 

1.00/401 

146.0 

1.47 (1.23-1.75) 






Joint 

2T 

1.20/937 

99.1 

1.02 (0.79-1.15) 

1.6 

49% 

38.9 (56.6) 

0.76 



mkcflow 

1.30/939 

40.0 

1.62 (1.53-1.83) 

4.7 

-0.78 

54.3 (37.6) 

0.68 

V1223 Sgr 

NuSTAR 

PL 

1.30/599 

233.5 

2.81 (2.30-3.41) 







mkcflow 

0.95/599 

222.5 

2.00 (1.67-2.39) 







Brems 

0.97/599 

119.5 

1.13 (0.91-1.38) 






Joint 

2T 

1.37/1141 

71.1 

0.62 (0.46-0.78) 

1.4 

48% 

22.2 (23.5) 

0.33 



mkcflow 

1.37/1143 

43.6 

0.66 (0.46-0.79) 

4.5 

-0.69 

34.3 (34.2) 

0.26 


^■866 text for full descriptions of the models. 

'^Difference in between fits with reflection amplitude let free to fit and fixed to 0.0 

‘^The column density for the partial covering absorber and the partial covering fraction (“2T”), or the 
maximum column and the power-law index for of the pwab absorber (“mkcflow”). 

‘^Temperature of the high temperature component (“2T”) or the highest temperature (“mkcflow”) in the 
best-fit model, with the best-fit value with reflection amplitude set to 0.0 in parentheses. 

® Overall abundance of accreting plasma, assumed to be responsible for both emission and reflection. 



